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Abstract 

We analyse supersymmetric models that show supersymmetry breaking in one and two dimensions using lattice meth- 
ods. Starting from supersymmetric quantum mechanics we explain the fundamental principles and problems that 
arise in putting supersymmetric models onto the lattice. We compare our lattice results (built upon the non-local 

C^ ■ SLAC derivative) with numerically exact results obtained within the Hamiltonian approach. A particular emphasis 

is put on the discussion of boundary conditions. We investigate the ground state structure, mass spectrum, effec- 
tive potential and Ward identities and conclude that lattice methods are suitable to derive the physical properties of 
supersymmetric quantum mechanics, even with broken supersymmetry. Based on this result we analyse the two di- 

r~^ ■ mensional N =1 Wess-Zumino model with spontaneous supersymmetry breaking. First we show that (in agreement 

with earlier analytical and numerical studies) the SLAC derivative is a sensible choice in the quenched model, which 
is nothing but the two dimensional 0"' model. Then, we present the very first computation of a renormalised critical 
coupling for the complete supersymmetric model. This calculation makes use of Binder cumulants and is supported 
by a direct comparison to Ward identity results, both in the continuum and infinite volume limit. The physical picture 
is completed by masses at two selected couplings, one in the supersymmetric phase and one in the supersymmetry 
broken phase. Signatures of the Goldstino in the fermionic correlator are clearly visible in the broken case. 
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> 

ff^ . 1. Introduction 

f I \ Symmetries are one of the guiding principles in contemporary theoretical physics. They led to the construction of 

('—^ ' the standard model of electroweak interaction \MM and strong interaction |4]. Based on the underlying symmetries 

bottom and top quark as well as the r neutrino have been predicted and the experimental discoveries, the last one 

more than two decades after the prediction |21, substantiated the success of the standard model to describe the physics 

on energy scales below 1 TeV. The extension of the standard model's symmetries by further ones is constrained by 

. . ■ the celebrated Coleman-Mandula theorem ||6(]. A way to circumvent this no-go theorem is given by extending the 

rS ' Poincare algebra with anti-commuting supersymmetry generators |7J that relate particles with integer spin to ones 

jrt ' with half-integer spin. 

The first field theoretical realisation of a renormalisable model with supersymmetry is the Wess-Zumino model 
Jsll in four spacetime dimensions with a field content of two (real) scalars and a Majorana fermion. Since then a 
variety of models with global supersymmetry have been constructed, for example supersymmetric gauge theories or 
supersymmetric sigma models and some of these possess several supersymmetries. For an introduction, see 191 41111 . 



Supersymmetric extensions of the well established standard model may solve or weaken several problems of the 



standard model, such as the hierarchy problem, the occurrence of dark matter and the strong CP problem 1112141411 . 

Supersymmetric models have certain theoretical advantages, such as less severe divergences in perturbation theory 
as compared to models without supersymmetry, and the supersymmetry algebra induces a vanishing ground state 
energy, as long as supersymmetry is unbroken. If the ground state is invariant under supersymmetry, mass degenerate 
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multiplets of bosonic and fermionic particles are predicted, and it known that in certain classes of supersymmetric 
theories a spontaneous breaking of supersymmetry is not possible 11511 . In experiments no such mass degeneracy has 
been detected so far, and the masses of bosonic and fermionic particles appear to be unrelated. At first si ght these 
results tell us that supersymmetry is not realised in nature. However, as first analysed by O'Raifeartaigh 11611 this 
non-degeneracy of masses is naturally expected if supersymmetry is dynamically broken. 

In a theory with dynamical supersymmetry breaking the ground state is not invariant under supersymmetry, and the 
ground state energy is lifted above zero jlTT . But the supersymmetry algebra is still present and this has implications 
for the physics in the dynamically broken sector. To date the Large Hadron Collider is fully operating and it is expected 
to measure remnants of supersymmetry in collision events within the upcoming years. Clearly, if supersymmetry plays 
any role in nature then it is mandatory to explore supersymmetric theories with methods that are applicable in the non- 
perturbative regime. 

Among those methods the lattice regularisation in combination with importance sampling based statistical 'Monte- 
Carlo' methods has been most successful over the last decades. Lattice methods often provide the only viable way 
to gain information about the non-perturbative sector of quantum field theories. Early simulations that aimed at an 
understanding of the pure SU{2) Yang-Mills theory JlSll built the basis for recent computations from first principles 
of the Hadron spectrum in full quantum chromodynamics 11911 . which is only possible due to increasing computing 
power and algorithmic improvements. As non-perturbative effects are automatically taken into account in lattice 
simulations, it is desirable to apply the lattice approach to supersymmetric theories as well. This has been the subject 
of a number of publications, see, e.g., i20l424ll and for recent progress in supersymmetric Yang-Mills theories 11251428 1 



and references therein. 

In all lattice regularised versions of field theories symmetries are of particular interest. If a symmetry of the contin- 
uum theory is not implemented in the lattice version it may happen that the symmetry is not restored in the continuum 
limit. E.g. for simulations of gauge theories it is important to implement the lattice version of the continuum gauge 
symmetry 12911 . But not every symmetry can be directly implemented in the lattice regularised theory. For instance. 



symmet 



the Nielsen-Ninomiya theorem II30I432I1 forbids the exact implementation of chirally symmetric fermions with a local 
fermion interaction and without introducing additional fermion flavours on the lattice. Nevertheless, it is possible to 
construct a (deformed) lattice version of the chiral symmetry, which is given by the Ginsparg-Wilson relation 1331 . so 
that a restoration of the continuum chiral symmetry is ensured in the continuum limit of the lattice action. 

For supersymmetry as extension of the Poincare algebra a fully realised supersymmetry algebra on the lattice must 
inevitably contain the generators of translations which would imply arbitrary translations to be part of the symmetry 
group of the lattice theory. By contrast, lattice regularised theories are only symmetric under translations by the lattice 
spacing. Therefore a complete realisation of the continuum supersymmetry algebra on the lattice is impossible and 
the full supersymmetry can only be realised as an accidental symmetry in the continuum limit of the lattice regularised 



theory. Technically, the reason for this can be traced back to the failure of the Leibniz rule on the lattice 113411 . 

It has been shown that even in supersymmetric quantum mechanics a naive discretisation does not lead to a super- 
symmetric continuum limit 13511 : generically, such a limit can at best be achieved by fine-tuning the bare coefficients 
of all supersymmetry-breakingcounterterms 13611 . This, however, requires much knowledge of the theory in advance. 
In some cases the relevant operators can be determined perturbatively, cf. 137!]. To circumvent the fine-tuning pro- 
cess several approaches are conceivable. Firstly a partial realisation of supersymmetry on the lattice is possible for 
theories with extended supersymmetry (for a review see 13811 ). Secondly recent developments aim at the construction 
of a Ginsparg-Wilson inspired relation for supersymmetric theories to obtain a lattice version of supersymmetry such 
that the continuum supersymmetry is broken in a controlled way i39il . Alternatively for scalar theories a deformed 
supersymmetry algebra on the lattice can be constructed by using a non-local product such that the theory is invariant 
under the full (deformed) lattice supersymmetry 1401 kill . 



Apart from an explicit supersymmetry breaking by the finite lattice spacing there exist further supersymmetry 
breaking effects that must be controlled in the analysis of supersymmetric theories. For example, at finite temperature 
Lorentz invariance and therefore supersymmetry as extension of the Poincare symmetry are brokenlJ In addition, 
for a finite spatial volume there may exist tunnelling processes between two formerly separate ground states such 



'This problem can be avoided by clioosing periodic boundary conditions also for tlie fermions. However, this is only possible for an unbroken 
supersymmetry. 



that the finite volume ground state energy is raised above zero. As it is inevitably to use finite lattices for numerical 
simulations these explicit supersymmetry breaking effects must be taken into account. 

We begin our investigations with discretised supersymmetric quantum mechanics with dynamically broken super- 
symmetry in Sec.|2] In this setting the basic concepts of supersymmetric theories are explained and reference results 
for certain observables are computed via the operator formalism, thus allowing to understand the physics behind su- 
persymmetry breaking on solid grounds. The corresponding lattice regularisation is based on a formulation that has 



been used in the unbroken supersymmetric quantum mechanics with great success II42I1 . We verify that even for quan- 
tum mechanical models with broken supersymmetry it is possible to obtain accurate results on the low lying energy 
spectrum from lattice simulations. 

The minimal setting for a field theory with supersymmetry breaking phase transition is given by the M = 1 Wess- 
Zumino model in 1 + 1 dimensions, which is analysed in Sec. [3] In the context of the quenched model a particular 
renormalised critical coupling for the Z2 symmetry breaking is shown to be independent of the chosen lattice regulator. 
The corresponding critical coupling in the full theory is determined and the relation between Z2 and supersymmetry 
breaking is worked out. 

2. Broken supersymmetric quantum mechanics 

An extensive analysis of quantum mechanical systems, such as the anharmonic oscillator, with lattice regularised path 



integrals has been performed almost three decades ago II43I1 . There has been renewed interest in quantum mechanical 



systems on the lattice in the context of supersymmetric quantum mechanics (SQM). In several works SQM has been 



used as a toy model to study the supersymmetry breaking induced by a naive lattice formulation 114411 and to explore 



lattice regularisations with partially 1351 WA |45L 14611 or fully 14011 conserved supersymmetries. It has been pointed 



out that a discretisation without any conserved supersymmetries may not be free of finite supersymmetry breaking 



renormalisation terms in the continuum limit 1I47I1 such that a careful treatment of supersymmetry restoration is needed. 
Most of the lattice studies of SQM so far have been carried out for the case of an unbroken supersymmetry and only 
few of them 1481 14911 consider the case of the dynamically broken supersymmetry. Here, the case of a SQM with 
dynamically broken supersymmetry is considered to explain the concepts and effects of supersymmetry breaking in 
a setting that allows for high precision measurements in the lattice theory and provides the possibility to compare to 
exactly calculable reference values from the operator formalism. 

2.1. Operator formalism 

SQM in one dimension is a generalisation of the supersymmetric harmonic oscillator|3 In analogy to supersymmetric 

field theories nilpotent supercharges Q and its adjoint Q) are introduced, 

with fermionic creation and annihilation operator '^^ and ^ and first order differential operators 

A^^ + P{<t>). At=-A+p(0) (2) 

containing the prepotential P{4>). In accordance to the field theory language we denote the position operator of the 
quantum mechanical system by 4>. The Hamiltonian is constructed via 

iiaan^CJ" ^^)^i(-^ + p'w + i*'.*iP'w).ff. (3) 

and acts on two-component state vectors |-0) = (|'0)b, |'0)f) where, for convenience, the first component is called 
'bosonic' and the second one 'fermionic'. The supersymmetry algebra is completed by the nilpotency of Q and Qt 
and the commutation with H, 

{Q,Q} = 0, {Qt,Q^} = 0, [Q,i/] = 0. (4) 



^An extended introduction to the operator formalism for supersymmetric quantum meclianics can be found in I50ll5lll 
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Figure 1 : Energy spectrum for the unbroken supersymmetric quantum mechanics discussed in 14211 . Q and Qt map between bosonic and fermionic 
sector. 



If P is a linear function of then A and A^ are the bosonic annihilation and creation operators of the (supersymmetric) 
harmonic oscillator. Accordingly the bosonic and fermionic Hamiltonian is given by 



H^ = \a^A = ^\^ + V^, i7p=l^At = -i-^ + yp, V^„^]-{P\cp)TP'm 



(5) 



Both Hamiltonians are by construction non-negative. The bosonic sector of a zero energy state is annihilated by A 
and a fermionic one is annihilated by A\ 



Hu 



^ A\0)b = 0, Hf\0)f = ^ At|o)p ^ 0. 



(6) 



The supersymmetry algebra implies a strict pairing of excited states, i.e. for every bosonic eigenstate lips) with energy 
E > there is a fermionic partner state 



I^f) 



1 



Q^I^b), I^b) = 



1 



/2E 



Qm 



(7) 



with identical norm and energy. 

The zero energy state(s) of the super Hamiltonian H can be given explicitly (in position space) as solutions of first 
order differential equations. If one of these functions is normalisable, then the supersymmetric ground state exists and 
supersymmetry is unbroken. Since the product of possible zero energy states (a;|OB) • (xjOp) is constant, there is at 
most one normalisable state with zero energy. The explicit form of solutions implies that for a polynomial prepotential 
P((/)) = J2n=o ^ri4'" with Cat 7^ supcrsymmctry is unbroken iff N is odd. In that case there is one normalisable 
zero energy state and the spectrum is similar to the one depicted in Fig.[Tl 

2.1.1. Supersymmetry breaking and the Witten index 

An existing and unbroken supersymmetry is defined by the existence of a normalisable grou nd state |0) which is 
annihilated by Q and Q^ which implies Hb\0)s = iJF|0)F = 0. Witten introduced an index JlSll to determine whether 
supersymmetry can be broken dynamically in supersymmetric field theories. In the present context it is given by the 
trace over all eigenstates of H, 

A = Tr(-l)^'', (8) 

where TVp = [ n i I is the fermion number operator that commutes with Hl^ Now, two alternatives exist: 



^As it stands, A is not well defined and requires a normalisation, e.g, A = lini^_>o Tr[e ^^ (—1)^^]. 
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Figure 2: Energy levels for the broken supersyimnetric quantum mechanics. Each level is doubly degenerate with one bosonic and one fermionic 
eigenstate. 



• For broken supersymmetry there is no normalisable zero energy state, 
energies and must be paired, which implies A = 0. 



All eigenstates of H have positive 



• For unbroken supersymmetry there are n^ bosonic and np fermionic ground states with zero energy. They 
contribute with n^ — np to the Witten index. All contributions from the excited states cancel, which gives 
A = riB — np- 

Therefore a non-vanishing Witten index implies an unbroken supersymmetry, but not necessarily vice versa. It is still 
possible that supersymmetry is unbroken while there are the same number of bosonic and fermionic zero energy states. 
For a one dimensional supersymmetric quantum mechanics at most one zero energy state is possible and A 7^ is 
equivalent to unbroken supersymmetry. 

2.1.2. Specifying the model 

The minimal modification of the supersymmetric harmonic oscillator with broken supersymmetry is given by the 

prepotential 



P{(t)) = mcf) + h(j)^ 



(9) 



with vanishing Witten index. Hence there is no normalisable ground state with zero energy. The spectrum is com- 
pletely degenerate and acting with the supercharges on one finite energy ground state will give the corresponding 
superpartner of this ground state 

This model depends on the dimensionful parameters m and h and in analogy to the supersymmetric harmonic 
oscillator m is used to set the scale. Therefore / = h/rn}-^ provides a scale independent dimensionless coupling. 
In a heat bath the dimensionless temperature is given by T = {m/3)^^ with (3 as dimensionful inverse temperature. 
Eventually coordinates are made dimensionless by setting $ = (p^/m. 

The energy spectrum and corresponding states can be computed directly by a discretisation of positio n sp ace 
and replacement of ^ by a discretised derivative. After the analysis of different possible discretisations in 15211 the 
numerically most stable choice is the SLAC derivative 15311 . which for periodic boundary conditions on a lattice with 



"^If not otherwise stated, "ground states" may also have a positive energy. 
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Figure 3: Probability density po{^) (siiaded areas) for bosonic and fermionic ground state and potentials Vg/p (lines) of the corresponding 
Hamiltonian at coupling / = 0.2. 



an odd number of points A^ and lattice spacing a is given by 



aid'^^^^y 



— (-Y 



sm(-K{x-y)/{Na)) 



N^ 



,2/aSLACN2 
Ixy 



-a\d^ 



N^ 3 

2tt^ / \x-v cos{TT{x-y)/{Na)) 

Ar2 I ) ^-^ 



sin-' {irix-y) /{Na)) 



: x = y 

■ x^y 

: x = y 

: x^y 



(10) 



The spectrum of the diagonahsed Hamihonian is depicted in Fig.|2] For weak couplings / < 0.1 there is an 
additional (approximate) degeneracy of the excited spectrum corresponding to the perturbed energy levels of two 
harmonic oscillators with energies Nm residing at the minima of the bosonic and fermionic potential V^^p (see Fig.[3]l. 

Interpretation as a physical system. The naming 'bosonic' and 'fermionic' sector may sound misleading because of 
the complete degeneracy of the spectrum. The system can be interpreted as a particle with spin 1/2 moving in an 
external potential that depends on the spin orientation. So 'bosonic' may refer to 'spin down' and 'fermionic' to 'spin 
up', respectively. Supersymmetry in this case is represented as degeneracy between an up and a down state. For 
the case of unbroken supersymmetry the ground state is unique and is invariant under application of supersymmetry 
although it is in a definite spin state given by the interaction potential. For the broken supersymmetry there are (in the 
present case) two different ground states none of which is energetically preferred (see Fig.|3]l. There is no interaction 
given by the Hamiltonian between bosonic and fermionic sector and one ground state will be preserved if no external 
interaction is applied (e.g. by interacting with a heat bath at finite temperature). Applying the supercharge will give 
the partner ground state and amounts to the symmetry between spin up and spin down state. Furthermore no linear 
combination of the two ground states is invariant under the supersymmetry. 

Physics flf T = 0. At vanishing temperature physics is given by ground state (vacuum) expectation values. Since 
supersymmetry is broken the system will stay in one of the degenerate ground states and expectation values are defined 
by this particular ground stateO Without loss of generality results are given for the bosonic (finite energy) ground state 
|0b) and the expectation value of an observable O is thus given by {0)q = (0b|O|0b)- 



'This is similar to the Z2 symmetry in the Ising chain. There, at any finite temperature the symmetry is restored. Only for T 
will take (and preserve) one of the possible "ground states". 



the system 



Observables can be computed from the diagonalised Hamiltonian to provide reference values for lattice compu- 
tations. This is done in the next section. On the lattice the primary focus lies on one- and two-point functions and 
the probability density of the coordinate $ given by po{^) = |('I'|0)| . The one-point function is then given by 

($)„ = /d$po(*)*. 

The bosonic two-point function (in the bosonic ground state) is defined through the Euclidean time evolution, 

($(t)$(0))o = (OB|$(t)$(0)|OB) = (OB|e*^$e-*^$|OB) = (ObI^c-'^^^^^'^IOb). (11) 

Equivalently the fermionic correlation function is computed by 

(*(t)vl/t(0))p = (OB|*e-*(^-^°)vl/t|OB). (12) 

In each case tm defines the dimensionless 'time'. 

The last quantity of interest is the effective potential which may be either defined by a Legendre transform of the 
Schwinger function^ or more directly at vanishing temperature by 

VM^o)= min {i'\H\^). (13) 

{^,|$|^)=$o 

Finite temperature physics. For any finite temperature there is a Boltzmann distribution with the same contribution 
of bosonic and fermionic partner states, including the lowest energy ground states. Again, with high precision calcu- 
lations of the low lying spectrum provided by the diagonaUsed Hamiltonian it is possible to compute the thermal field 
distribution and expectation values 

E E E 

where the sums run over all bosonic and fermionic states. 

2.2. Lattice regularised path integral 

With the methods given in the previous section it is possible to obtain exact results against which the path integral 
based calculations can be compared. Therefore the accuracy of the lattice simulation can be determined even when 
supersymmetry is broken. The corresponding Euclidean path integral is given by 

Z= fv(j)Vi'V'iije-'^^'^-''f'-'f'\ (15) 

with Euclidean action 

^ = ydT Q(50)2 + ip2(^)^^(5^p,(^)) A (jg) 

Expectation values are computed via 

{A)=Z-^ I V(l)VipVi)A[(f),iP,i!]e-^^'l'''''''l'\ (17) 

Supersymmetry appears as a symmetry of the action, where one transformation is given by 

5^^U = H', <5(i)i/i = -£((/) + P((/))), (5(iV' = (18) 

and a variation of the action gives j'^-'S' = /dr [d^ePtp)] = 0. In the same way the action allows for a second 
supersymmetry transformation 



The Schwinger function is naturally defined in a path integral formulation. 

7 



For the above supersymmetries to hold it is necessary that the fields vanish at infinity or that they are periodic in the 
EucUdean time. But for a thermal path integral at inverse temperature /3 the fields obey the boundary conditions 

0(O) = 0(/3), V(0) = -^(/3), V^(0) = -V;(/?), (20) 

i.e. the fermionic field is antiperiodic in time. Since the fields need not vanish anymore the variation of the action 
then reads 

<5(i)^ = [ePV]f.o = -2 [^"^V'].=o (21) 



which can be non-vanishing so that supersymmetry is broken by the finite temperature. In 11421 14411 for an unbroken 
supersymmetric quantum mechanics periodic boundary conditions have been used to avoid such an explicit breaking. 
For temperature going to zero a change in boundary conditions is equivalent to an insertion of (—1)^'' into the path 
integral, 



Here, the periodic path integral is vanishing due to A = for a broken supersymmetry. Thus, for a theory allowing 
supersymmetry breaking, periodic (supersymmetry preserving) boundary conditions cause a severe sign problem. 
This does not completely rule out the choice of these boundary conditions, as will be discussed on the case of the 
two dimensional Af ^ 1 Wess-Zumino model in Sec. 13.21 but puts constraints on the range of applicabilityO To have 
a well defined (non-vanishing) path integral antiperiodic (thermal) boundary conditions for the fermionic fields are 
imposed. 

For a construction of a lattice model the choice of the lattice regularised derivative is crucial^ The canonical 
choice for scalar theories would be the forward (or equivalently backward) derivative. For derivatives appearing in 
the fermionic action a popular choice is given by Wilson's prescription 12911. Nevertheless, these simple discreti- 



sation rules are not applicable to supersymmetric theories as analysed in 11421 l44l 14711 for the case of an unbroken 



supersymmetric quantum mechanics. These results show the need for a more careful treatment of the discretisa- 
tion of supersymmetric theories. In the comparative study of six different discretisations li42il the one based on the 
SLAC derivative provides results close to the continuum limit even at finite lattice spacing. For an odd number of 
lattice points with periodic boundary conditions the matrix representation is already given in Eq. ( fTOl i. Antiperiodic 
boundary conditions (necessary for fermionic fields) are best realised on an even lattice with A^ points, 

«(^ )^y=')jL(\{=^-y)/a 1 -t^h' '■^■^^ 

{N^ I sin(TT(x-y)/{Na)) ■ -^ r y 

while the squared SLAC derivative for an even number of lattice points and periodic boundary conditions (as needed 
for the bosonic fields) reads 

2/aSLAC\2 J #^" 3' ~ : X = y 




a^(5^^^x = <!^;^.^r.-,)/a 1 .,[■ (24) 

sin^{-n{x-y)/{Na)) ■ -^ r U 

Although it was analysed 15411 th at this prescription will lead to a non-covariant and non-local continuum limit in 
lattice QED it can be proven 11421 15511 that for scalar theories in one or two dimensions with Yukawa interactions a 
local renormalisable continuum limit is reached. For that reason the SLAC derivative is used here to regularise the 
supersymmetric quantum mechanics on the lattice with corresponding action 

^ = - E U^^d'^^^^f^y^y + 1 E P^^^f + E ^- (^'^^^ + P\^^)5^y)^y (25) 

x^y X x,y 

on lattices with an even number of sites, where field (p and derivative d^^^'^ are dimensionless and arise from a 
rescaling of the dimensionful quantities with the lattice spacinglj 



'For a phase with unbroken supersymmetry in models with A = only one specific ground state belongs to the physical spectrum and periodic 
boundary conditions may be imposed. 

**In contrast to the operator formalism where the field space is discretised, the lattice path integral is based on a discretisation in the Euclidean 



time. 

9' 



The fermionic fields ip, ip are already dimensionless and do not need to be rescaled. 
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Figure 4: Sign of the fermion determinant measured at fixed / = 1 (left panel, 10*" configurations per data point) and at fixed mP = 36 and 
N = 50 lattice points (right panel, 5 ■ 10^ configurations per data point). 



Access to the non-perturbative sector of the lattice model is gained from Monte-Carlo simulations which have be- 
come a powerful tool due to increasing computer power and algorithmic improvements that allow for the inclusion of 
dynamical fermions in simulations. These statistical methods are based on importance sampling and the interpretation 
of the lattice regularised path integral 

Z= /"D0DV'DV;e-^['^''^''''l (26) 

as probability distribution. To construct the probability density the action is then split into a bosonic and fermionic 
part according to 

S[^, V, ^] = 5b [0] + ^ i^,M,y[^]ijy. (27) 

Applying the rules of Grassmann integration the fermionic part of the path integral can be integrated out and yields 



Z= B4> detM[(^]e"^»l'^l. 
In this way (bosonic) expectation values are computed by 



In a Monte-Carlo simulation the lattice regularised fields cj) are generated according to the distribution 



pM 



^-SB[</']+ln dctA/[0] 



(28) 



(29) 



(30) 



After a number of iVMc samples one obtains a time series (j)'^''\k ~ 1 , . . . , ^Vmc, and expectation values are evaluated 
using 



(O) 



Nmc 
fe=l 



(31) 



This expression is only exact iff det M > 0. If det M is negative the sign has to be taken into account by reweighing. 
However, the emphasis shall lie on the physical questions and for further simulation details the reader is referred to 



the rich literature on Monte-Carlo methods, e.g. II56M61I. 
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m 



2.2.1. Sign of the fermion determinant 

For periodic fermionic boundary conditions Zp ex A will vanish in the continuum and a severe sign problem is 

expected to arise in reweighed expectation values. For thermal boundary conditions it is a priori unknown if there are 



p(*) 
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Figure 5: The distribution of the averaged field # for mf} = 4 at coupling / = 0.2 with respect to the sign of the determinant for periodic (left 
panel, A'^ = 101) and antiperiodic (right panel, A^ = 100) fermionic boundary conditions obtained from 10^ configurations. 

configurations with dot M < and if there is any dependence on lattice spacing, temperature, or coupHng. For that 
reason (sign det M) has been measured in the sign quenched ensemble with the distribution given by Eq. ( l30l l for 
different parameter sets (see Fig. |4]i. These results imply a complete absence of the sign problem in the continuum 
limit for every coupling and temperature. The sign problem only exists for large couplings / at fixed lattice spacing 
and temperature. 

2.2.2. Ground state structure 

With thermal as well as supersymmetry preserving (periodic) boundary conditions for low temperature T = 0.25 the 
ground state structure is analysed!^ Simulations at / = 0.2 are performed and the distribution of the lattice averaged 
field $ = A^~^ ^^ $2; is analysed with respect to the sign of dot A/ (see Fig. |5]l. Configurations with $ > —1/(2/) 
are unaffected by a change of boundary conditions whereas the sign of det M changes for $ < —1/(2/). This 
behaviour can be seen explicitly on the level of the discretised action. For the chosen prepotential bosonic and 
fermionic ground state are related by a Z2 symmetry ^^ — > —^x ~ 1//- S'b is invariant under the symmetry 
operation whereas the effect on the fermionic contribution depends on the derivative used. The SLAC derivative has 
an antisymmetric matrix representation, d^]:^'^ — -9^^^*^. P'{4') enters on the diagonal of the fermion matrix M. 
Applying the Z2 symmetry gives P'{<t>) -^ —P'{4)) and changes the sign of the diagonal elements of the fermion 
matrix. Altogether, the symmetry operation changes M{(p) — > — A/^(0). For antiperiodic (periodic) fermions the 
fermion matrix size will be even (odd, respectively) and the determinant will keep the modulus but changes its sign 
for periodic boundary conditions while for antiperiodic fermions the sign is preserved. Therefore periodic SLAC 
fermions imply Zp = exactly. The boundary condition dependence of the distribution coincides with introducing 
(— l)^"" into the path integral for periodic boundary conditions and configurations with $ > —1/(2/) correspond to 
the bosonic ground state whereas the other ones correspond to the fermionic ground state, respectively. This is in 
accordance with results from the operator formalism shown in Fig.|3] 

2.2.3. Thermal field distribution 

At finite temperature the single site distribution pt{^) of Eq. (fT4l i is computed on a lattice with N = 100 points 
at fixed coupling / = 0.2 in the temperature range m/3 e [0.2,4] with 10^ configurations. Even for this coarse 
lattice the results match almost perfectly with the reference values from the diagonalised Hamiltonian (see Fig. |6] 
left panel). At low temperature m/3 = 16 it is possible to extract the probability distribution in the bosonic ground 
state by considering only configurations with $ > —1/(2/). Even at finite (but small) temperature the so-extracted 
probability distribution follows the exact result at T = 0, see Fig. |6] (right panel). In consequence, for T ^ a 
thermal mixture of bosonic and fermionic ground state is found whereas at exactly vanishing temperature the system 
can be triggered to stay in one ground state. 



'"Periodic boundary conditions have only been used for this particular analysis of the ground state structure. 
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Figure 6: Probability distribution p{^) for different temperatures at coupling / = 0.2 on a A'^ = 100 lattice. Lines depict the exact results, points 
arise from the lattice calculations. For compaiison the distribution at T — > is also drawn. Left panel: Thermal distiibution. Right panel: Exact 
distribution for the bosonic ground state and distribution measured on configurations with <!> > 



1 

'2f 



2.2.4. Effective potential 

The effective potential as introduced in Sec. 12.1.21 is not directly accessible in lattice simulations. 

quantity is given by the constraint effective potential 16211 



f/(^o) 



1 

1 



In ( / D0 dot M[4)] c-^''['^I 5(l> - $0 



A better suited 



(32) 



which can be straightforwardly computed on the latticel^ It has been proven for the case of bosonic scalar theories that 
the constraint effective potential coincides in the limit of infinite spacetime volume (or here equivalently, vanishing 
temperature) with the effective potential, lim,„^_i,oo U{^) = lim,„^_>oo K:ff (*&)■ 

We determined the constraint effective potential for various inverse temperatures m/3 at fixed coupling / = 0.2 
for N = 300 lattice points. The result is depicted in Fig. [7] together with the effective potential for the full theory and 
for the bosonic/fermionic sector To avoid the ambiguities of the unknown additive constant, the minimum of each 
potential is normalised to OPi 

It is often stated that the conventional effective potential for a quantum mechanical system is strictly convex. This 
is true for systems with Hilbert space L^(M''). For the supersymmetric system the Hilbert space is i^(M) (g) C^ and 
all linear combination of the bosonic and fermionic ground states have the same energy, which implies a flat region in 
the effective potential as defined in Eq. ( fT3l l. 

To compute the value of the constraint effective potential even around the peak $ = —1/(2/), reweighting 
methods ll63ll similar to the multicanonical ensemble l64ll have been applied. Simulations were performed with a 
modified action S'g = Sb ~ W{^), with VF($) suitably chosen by iterative refinement to approximate /3J7($) 
between the two minima and W^($) constant in the outer regions. Of course, the configurations must be (re-)weighted 
after simulation with a factor exp(— VF($))o 

It is apparent that the positions of the minima of L/($) correspond to the bosonic and fermionic ground state. How- 
ever, the constraint effective potential shows no tendency to flatten out towards the conventional effective potential. It 
rather tends towards 

[/(!>)= lim minfl/eff,B($),Kff,F($)), (33) 



lim 

m/3— fo- 



lim 



where T4ff_B/F denotes the effective potential of the bosonic (fermionic, respectively) sector of the Hamiltonian (|3]l. 

This result can be explained by the behaviour of fermions on the lattice under a change of boundary conditions. 
The ratio of the fermion determinant between periodic and antiperiodic boundary conditions can be computed in the 



' ' Only an additive normalisation constant of the constraint effective potential is left undetermined. 

'^This amounts to an offset for the effective potential of the full theory that is given by the ground state energy, AVJ^ff = Eg = 3.68 • 10~^m. 

'^^In a conventional Monte-Carlo simulation without reweighting the configurations around the peak will be suppressed at m/3 = 49 by more 
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Figure 7: Effective potential for the full model ("exact fuU"), bosonic sector ("exact bosonic"), and fermionic sector ("exact fermionic") of the 
Hamiltonian at T = and constraint effective potential from lattice simulations of the full model at various m/3 for coupHng / = 0.2 and 
A'^ = 300 lattice points, measured with up to 10^ configurations. 



continuum limit Il65n . and by using zeta function regularisation 116611 one gets rid of the free determinant as normalisa- 
tion constant. The continuum result 



R U] ^ detAfp[</.] ^ sinh[l/pV^(0(r))dr] 

^^^^ det Map[0] cOsh[i /„^ P'(0(T))dT] 



(34) 



will be either +1 or —1 in the infinite volume limit, depending on the value of $, 

1 f~i ' 

2-' 



lim Rp{(j)) = 



$ > 
$ < 



(35) 



which is equivalent to (—1)^''. I.e by means of the lattice path integral with fixed (thermal) boundary conditions it 
is only possible to assess either the zero fermion sector or the one fermion sector. The interpolating states between 
both sectors that are responsible for the flattening of the effective potential are thus not accessible by the (lattice) path 
integral at any finite temperature. 

2.2.5. Two-point functions and spectrum 

Lattice based path integral methods provide a non-perturbative way to gain information about the spectrum of the 

theory. Via the long distance behaviour of correlators C{t) >■ exp(— mphysi) it is possible to extract the physical 

'pole mass' TTiphys which is given by the imaginary part of the pole of the propagator G{p) = {FC){p), the Fourier 
transform of the correlator, and describes the energy difference between the ground state and first excited state of the 
theory. For that reason connected correlation functions in the thermal ensemble have been computed for bosonic and 
fermionic fields (see Fig. [8]) with 



CB(i) = ($(t)$(0))-($)' and CF(t) = (V(t)^(O)) 



(36) 



The correlators take non-vanishing constant values for large distances in a region where the exponential falloff drops 
below the visibility scalel^ Correlators computed from the lattice regularised theory fit nicely to the ones computed 



'"^This is not directly related to any unconnected part. Here only connected correlators are considered. 
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Figure 8: The two-point function (left panel: bosonic, right panel: fermionic) for the thermal ensemble given by m/3 = 9 and coupling / = 1 
computed by the diagonahsed Hamiltonian and on a lattices with N S {100, 200, 400} points. For the bosonic case the data points coincide with 
the reference values. A statistics of up to 10^ configurations has been used. 



Csit) 





Figure 9: Bosonic (left panel) and fermionic (right panel) correlator in the ensemble projected to one ground state as obtained from the diagonahsed 
Hamiltonian and from lattice simulations with A^ = 400 sites at coupling / = 1 and statistics of up to 3 ■ lO'^ configurations. 

by the diagonalised Hamiltonian. Fluctuation are still visible around the continuum values with the size of fluctuations 
vanishing for smaller lattice spacings. Further it is possible from the fermionic correlator at small (but non-vanishing) 
temperature to compute the overlap of bosonic and fermionic ground state I (Of|^^ |0b) I = 0.41174 by the approxi- 
mation CF{m(3/2) = 0.419(7) for m/S = 16 and N = 400. 

With nearly vanishing temperature the system will mainly reside in the ground states, and with the results of 
Sec. 12.2.21 it is possible to compute the correlation function at T = by projecting to one of the ground states 
(see Fig. |9]). In that case the bosonic correlator shows no constant part and the exponential behaviour completely 
coincides with the one resulting from the first excited (bosonic) state. From the lattice values at A^ = 400 a mass of 
'TT'phys = 1.5047(15) m is extracted through an exponential fit in the range tm G [1, 3] from a simulation at 77if3 = 25. 
This is very close to the exact value Ei ~ Eq ~ 1.5046 m. Of course this method is only applicable for t < /3/2 and 
there will be large deviations for t close to /3/2. 

All these results clearly demonstrate that it is possible to extract correlators at finite and zero temperature from a 
lattice discretisation in complete agreement with exact results in the continuum limit. The degenerate ground states 
are visible as a constant part in the fermionic correlator and a projection to one ground state allows to extract the 
energy difference between the ground state and first excited state (see Fig.fTOll. 
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Figure 10: Excitations visible as t — > oo behaviour of the correlators evaluated on configurations projected to the bosonic ground state. 



T 



T~ 



exact 



extrapolated 



1.00 
1.44 
4.00 



0.4392039 
0.4027741 
0.3730929 



0.43902(18) 
0.40298(20) 
0.37327(18) 



Table 1: Linearly extrapolated lattice values of the Ward identity (P) /\/m in comparison to the exactly calculable ones from the operator 
formaHsm for three different temperatures T at coupling / = 1. 

2.2.6. Ward identity 

For field transformations $' = <1> + (5$ that do not change the path integral measura'^l (V^ ~ P$') Ward identities 

naturally arise on the level of observables as 

(O) =Z-^fv<^'0[<^']e-'^^'^'^ =Z-i/x>$(0[$]+(50[$])(l-(5S'[$])e-^[*l = {0 + 50-05S), (37) 

implying (SO) = {O SS) . If further the action and ground state are invariant under the transformation given by 6 then 
((50) vanishes for every observable O. 

For unbroken supersymmetric theories Ward identities are used to test the supersymmetry restoration in the con- 
tinuum limit by analysing the continuum limit of (SO) for a given observables O. If supersymmetry is broken then 
Ward identities will not be fulfilled in the continuum limit, {50) ^ 0. On the lattice supersymmetry will be further 
broken explicitly by a finite lattice spacing and by finite temperature. 

A simple Ward identity is provided by 



di 



(<5(i)i.) = -£ 



di 



+ P(<^)) = -£-( / diF(0) ) = -/3e-(P(0)) , 



(38) 



and a dimensionless identity is given by (P) j \fm = iff the ground state is invariant under the supersymmetry. The 
impact of the explicit supersymmetry breaking on the continuum result at vanishing temperature is depicted in Fig. [TT] 
for coupling / = 1. The finite a effects are small for the considered lattice spacings and a linear extrapolation to a = 
provides results in full agreement with the continuum results (see Tab. [T]i. Finite temperature effects are negligible 
for T < 0.2. The numerically exact value from the operator formalism is (P) /\/rn = 0.37251 and corresponds to 
the non-supersymmetric ground states. This quantity is invariant under the Z2 symmetry $ ^ $ — 1// and is not 
sensitive to the specific ground state chosen at T = 0. 



"Here only the anomaly free case is considered. If the path integral measure is changed under the transformation additional contributions must 
be taken into account. 
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Figure 1 1 : Prepotential for coupling / = 1. Left panel: Continuum limit of the prepotential (P) I \/ra which serves as simplest Ward identity for 
the supersymmetric quantum mechanics. (The dashed lines denote the continuum values while the shaded area gives the en'or bound of the linear 
extrapolation to a = 0.) Right panel: Temperature dependence of of (P) / \/rn computed from the diagonalised Hamiltonian (data from the left 
panel is marked by points). 



3. A/" = 1 Wess-Zumino model 



Supersymmetry as incorporated in the minimal supersymmetric standard model 111211 or extensions thereof II67I1 can 
only be a fundamental symmetry of nature if it is spontaneously broken on experimentally accessible energy scales 
with a phase transition at a much higher energy. Guided by this observation there is a need to study supersymmetry 
breaking phase transitions with non-perturbative tools. To explore the possibilities provided by lattice regularisations a 
minimal model is chosen, namely the N = (1,1) Wess-Zumino model in two dimensions [68]. It is obtained from the 
N = (2, 2) Wess-Zumino model, which provides a dimensionally reduced version of the matter sector of the minimal 
supersymmetric standard model, by constraining the fields of the M = (2, 2) model to be real JA'A kg, l60l 169^7 111 . 
This amounts to turning complex scalars into real ones and replacing Dirac fermions with Majorana fermions, leading 
to a minimal field content with only one bosonic and fermionic degree of freedoml^ For models with short range 
interactions at least two dimensions are necessary to observe a phase transition. In particular supersymmetric quantum 
mechanics is not sufficient to model a phase transition since the systems reside in one specific phase depending on 
the highest power of the superpotential (cf. Sec. 12.11 ). Since the seminal work by Witten 11511 it is known that the 
index Tr(— 1)^'' can vanish for specific choices of the prepotential and supersymmetry may be broken spontaneously 
depending on the couplings of the prepotential for the M = \ Wess-Zumino model. 
The off-shell continuum formulation is given by the action 



S ■ 



<^^x-{{d,^f 



- 4}{$ + P'(0))V + 2FP{(j}) - i^2 



(39) 



where F and denote real scalar fields and ?/; is a (real) Majorana spinor with two spinorial components. This 
formulation is invariant under the supersymmetry transformations 

(40) 



(41) 



(50 = eiA, 5%l; = (^0 - F)e, 5F = -e^V- 
By eliminating the auxiliary field F = P{(p) one arrives at the on-shell action 

s= U^x\ {{d^^f + ii$ + p'{<p))4> + P(0)^) . 



Here, the prepotential P is chosen to be 



P(0) 



A^o 



V2A V 2^ ' 



(42) 



'^From a practical point of view the absence of gauge fields has the advantage that derivatives can be applied in momentum space to speed up 
Monte-Carlo simulations considerably. 
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Figure 12: The classical bosonic potential of the Wess-Zumino model given in Eq. )43t . 



SO that the classical potential (Fig.[T2]l for the scalar part corresponds to a (p'^ theory with Z2 symmetry (0 



-'^), 



PW' 



Mo ,2 



const. 



(43) 



while the fermions acquire a Yukawa interaction. An analysis of the Witten index 11511 reveals one bosonic and 
one fermionic ground state that imply Tr(— l)^"" ~ 0. This ground state structure then allows for a spontaneous 
supersymmetry breaking. It is expected that for fixed A > and /Iq <C the system cannot tunnel between the two 
ground states so that supersymmetry is unbroken. On the other hand, for jIq > both ground state energies are lifted 
above zero and supersymmetry is broken. In consequence there will be a supersymmetry breaking phase transition 
for some /i§ < and this section is devoted to the study of this transition. 

There have already been several studies aiming at analysing and understanding supersymmetry breaking in this 
model. By a strong coupling expansion 117211 the occurrence of a supersymmetry breaking phase transition was pre- 
dicted. Calculations of the ground state energy with Monte-Carlo methods ll73ll confirm this expectation. Afterwards 
lower bounds on the ground state energy have been analysed 1174 17511 to obtain a phase diagram of supersymmetry 
breaking by working in the Hamiltonian formalism and making a numerical analysis with Green's function Monte- 
Carlo methods. However, the obtained critical lattice couplings so far are only unrenormalised couplings correspond- 
ing to one specific lattice spacing. Recently an analysis based on exact renormalisation group methods ll76ll has been 
performed ll77ll where no supersymmetry breaking is introduced during the renormalisation group flow. Again, the 
phase diagram has been obtained and supersymmetry breaking is found to coincide with a restoration of the Z2 sym- 
metry in a second order phase transition. A supersymmetric (massive) phase is found for weakly coupled systems 
as well as a phase of broken supersymmetry with goldstinos, the massless Goldstone fermions 117 811 of the broken 
supersymmetrjl^ and bosons whose mass vanishes with growing renormalisation group scalej^ But still the critical 
coupling depends on the chosen regulator which prevents a direct comparison of numerical values. 

In general a naive discretisation of a lattice action has to face the problem of broken supersymmetry for finite 
lattice spacing with the need to fine-tune lattice couplings to reach a supersymmetric continuum limit. However, 
this model has the advantage that the necessary counterterms have been analysed in lattice perturbation theory 137 1 
and a lattice prescription is given that ensures the supersymmetric continuum limit (perturbatively). Simulations of 
the given discretisation (based on the Wilson derivative) have already been performed JSOll and a tunnelling between 
the possible ground states is found to coincide with the onset of supersymmetry breaking and the appearance of a 



"This breaking is not forbidden by the Mermin- Wagner theorem 17911 that only applies to bosonic symmetry generators. 
"*Also the corresponding critical exponents have been determined in 17711 . 
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Figure 13: The only divergent Feynman diagram for the bosonic <j>^ model in the Z2 symmetric phase. 

goldstino. Despite all these successes several open issues remain: The breaking of a Z2 symmetry (which is correlated 
with the restoration of supersymmetry) has so far only been analysed with a Z2 breaking actionl^ Further the given 
critical coupling is still regulator dependent and not directly comparable to other methods. For that reason the aim 
of this section is the non-perturbative determination from first principles of a renormalised critical coupling in the 
continuum limit. 

3.1. Quenched model 

Although our primary focus lies on the supersymmetric model it is useful to exemplify the definition of a renormalised 
critical coupling in the setting of the quenched model where fermionic contributions are neglected. This becomes 
especially important because a non-standard discretisation based on the SLAC derivative is used. 

With the prepotential of Eq. ( |42] | the quenched model is nothing but the two dimensional 0^ model with action 

This model is (classically) invariant under a discrete Z2 symmetry (0 -^ —(f>) which can be broken dynamically in 
the full quantum theory 18 ill . The symmetric phase is characterised by (0) = 0, whereas in the broken phase (in the 
thermodynamic limit) (0) 7^ 0. 

In contrast to the full Af ^ 2 Wess-Zumino model the (f)"^ model (as well as the full J\f = 1 Wess-Zumino model) 
is not finite and there is need for a renormalisation of couplings. In the Z2 symmetric phase the only divergence arises 
from the 'leaf diagram (see Fig.fTSll and the model can be made finite with a mass renormalisation. 

The (one-loop) relation between the inverse propagators is given in terms of renormalised couplings by 

G-\p)=p' + fi' + j:{p^), J:{p') = 3XA^2-S^i', A^2= [ '^P ^] ^ (46) 

J [i.-Ky p-^ + /i^ 

and the only (logarithmically) divergent contribution is A^2 . Hence all ultraviolet divergences can be removed by a 
renormalisation of the mass terno 

(5/^2 = 3AA^,2 ^ fil = 11^ -Z\A^2. (47) 



This renormalisation prescription will also be sufficient in the Z2 broken phase 118 lH . We follow f820 and use the 
dimensionless renormalised coupling / = X/ n^ that distinguishes between the symmetric and broken phase. By 
contrast, (naively) possible definitions of a renormalised coupling that are based on the pole mass or the propagator at 
vanishing momentum are tiot sufficient to distinguish between both phases. Note that this renormalisation prescription 
can only be applied in a given scheme and continuum results will then follow by removing the ultraviolet regulatorr'l 



"The Wilson term for the fermionic part of the action will break the Z2 symmetry as analysed in I6OI1 . 

^"Equivalently this renormalisation can be obtained by a normal ordering of the interaction part in the symmetric phase, 



s, = |d^xl(^i^,4>f + ^^^4>' + ^■.^^■.,y 



where : . : ^ denotes normal ordering with respect to the mass /j. 
^' There is no need to introduce an infrared regulator. 
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Figure 14: Extrapolation of the critical couplings of the bosonic (/>* model to the continuum hmit. The shaded region indicates error bands for the 
extrapolation using the functional form 152) . Data points are taken from I83n with the canonical discretisation while confidence bands are computed 
here. 



3.1.1. The 111 phase transition on the lattice 

The computation of the critical coupling in the bosonic 0* model has a long history where several methods and 
approximations (e.g. based on the Gaussian effective potential or light-cone quantisation) have contributed. For a 
recent overview containing the most precise Monte-Carlo results see 18311 . Since these accurate results serve as 
reference values in the present work we briefly recall the most important points. The lattice regularisation in 182118311 
is based on the "canonical" discretisation utilising the forward derivative and a renormalised critical coupling is 
computed with high precision. The canonical model is given by 



Sv, 



eMe( 



*x)' + Ao'^' + 2' 



(48) 



with dimensionless lattice parameters A = \a} and /Jq = [i^o^ ■ Since A acquires no renonnalisation it is used to 
set the scale. Equivalently A determines the lattice spacing with A — > in the continuum limit. The (dimensionless) 
renormalised coupling is again given by / = X/jl^ and the corresponding jj.^ can be computed via /Iq = /t^ — 3\Aj^2, 
where A^2 that enters the normal ordering is given in the infinite volume limit by the lattice propagator for the forward 
derivative, 

n n _. 

Aa2 = lim n"^ y^ y^ _ /49) 

^ "^oo ^^ ^^ /i2 ^ 4 gin^ (Trfci /n) + 4 sin^ (Trfcz/n) 

Because A^2 diverges only logarithmically for a ^- it follows that fi^, fi^ -^ in the continuum limit at fixed 
/. One can show that at every fixed A a second order phase transition appears. However, the continuum physics at 
these phase transitions corresponds to an infinite continuum coupling A and is therefore only indirectly relevant for 
the continuum 0^ model under consideration. 

At finite lattice spacing (given by fixed A) there will be a Z2 breaking phase transition and a critical fi^ can be 
extracted. The renormalised critical coupling in the continuum hmit is then determined via 
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Figure 15: Binder cumulants at A = 0.15 for different lattice sizes. The coarse view is given in the left panel while a close-up view of the critical 



region is shown in the right panel. From the intersection point for the largest lattices a critical coupling /c 
extracted. The shaded region gives the error band for the infinite volume extrapolation. 
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0.01325(5) 



11.321(43) is 



The phase transition itself for finite A can be determined from the Binder cumulant 



U = 1 



All. 



with 



N-'Y.'t'-^ 



(51) 



which becomes independent of the lattice volume at the second order phase transition point 18411 r^l The critical 
coupling has been determined from lattices up to a size of 1200^ and was found to be affected by linear and logarithmic 
corrections in the lattice spacing. An extrapolation based on the published values in ll83ll for A e [0.01, 1] to the 
continuum using a functional form 

/c(A)«/c + aA + 6AlnA (52) 

reveals a renormalised critical coupling in the continuum of /c = 10.81(7) (see Fig.[l4]i. 



3.1.2. Regulator independence of the renormalised critical coupling 

The results of Sec.|2]and results in ll40ll42, 46, 60, 7111 imply that a discretisation based on the SLAC derivative gives 
results close to the continuum limit and will not break the Z2 symmetry of the full supersymmetric model. For that 
reason the SLAC derivative will also be used to simulate the N = 1 Wess-Zumino model. In contrast to other low 
dimensional scalar supersymmetric models it is now necessary to cope with a logarithmic mass renormalisation and 
the renormalised lattice coupling depends on the chosen renormalisation procedure. 

To justify the applicability of the SLAC derivative also for the present case simulations based on the lattice action 



^ = eUE(^^^'<^)'+a^^' + ^ 



(53) 



have been performed where the mass parameter is still given by /ig = /i^ — 3AA^2. Only now, A^2 is determined 
from the propagator based on the SLAC derivative, 



^^2 = lim {2ny 



^ E E 

fci— — n+l A;2— — n+ 



^ /i2 + (7rfci/n)2 + {■Kk2lnY 



(54) 



Similar to the case of the naive discretisation the crossing of the Binder cumulant U for different lattice volumes 
at fixed A determines the critical /t^ and therefore the critical coupUng /c. This procedure is exemplified for A = 0.15 



■^^Strictly speaking there is still a slight volume dependence such that the large volume extrapolation of the intersection points for different lattice 
volumes corresponds to the phase transition. 
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A 


/c 


0.02 


11.035(76) 


0.05 


11.112(74) 


0.10 


11.268(57) 


0.15 


11.321(43) 


0.20 


11.386(42) 


0.25 


11.429(29) 



Table 2: Renormalised critical couplings for the <f>'^ model as determined from lattice sizes up to 256^ with the SLAC derivative. 

in Fig. [TSJ where lattice sizes up to 256^ were used. As an outcome of these calculations critical couplings have been 
determined for six different A (see Tab.|2]i. Again, an extrapolation to the continuum limit A ^^ has been performed 
according to Eq. ( |52] | (see Fig.fTSIl and gives the renormalised critical coupling in the continuum of /^ = 10.92(13) 
which is in complete agreement with the reference value /c = 10.81(7) of 18311 . This proves that the lattice t/)^ 
model with non-local SLAC derivative possesses the correct continuum limit, as expected from the analytic studies in 



1 42l |55|] . The definition of a renormalised critical continuum coupling is independent of the chosen lattice regulator. 



3.2. Full dynamical model 

By inclusion of dynamical fermions the model is now given in the continuum by the action W\\ . Irrespective of the 
chosen prepotential P((/)) the action is invariant under one supersymmetry. The dynamical breaking of this supersym- 
metry for the prepotential defined in Eq. (l42l l will be analysed in the following. 

3.2.1. Renormalised lattice parameters 

Using a discretisation based on the SLAC derivative the lattice action is given by a direct discretisation of the corre- 
sponding continuum action, 

S = T.\ (E(^''^'^)' + Ag^l + \^l + ll V^I C(C^+V2l^^ V^,) . (55) 

^ =A/[0] 



2 

X V 



In a Majorana representation with 7 matrices and charge conjugation matrix 



-1; _»,V .M\ -X c-{\ -, 



is the fermion matrix M real and antisymmetric. 



In 13711 a lattice model of the M ~\ Wess-Zumino model built upon the symmetric derivative has been analysed. 
To avoid the species doubling problem a Wilson mass term has been added to the prepotential. Similar discretisations 
were used for the M =2 Wess-Zumino model in 1421 l46l rod ItoII? Ill . With properly renormalised coupling parameter 



(as described below) is was shown that the supersymmetric continuum limit is reached. Although this restoration 
was only analysed for a discretisation based on Wilson fermions the whole line of argument is directly applicable for 
the SLAC derivative. This follows from the fact that the SLAC derivative for two dimensional models with Yukawa 
interaction needs no non-local or non-covariant counterterms to construct a local continuum limit and the lattice 
degree of divergence of Feynman diagrams coincides with the divergence of the corresponding continuum diagrams 

iSSllil. 



The analysis of divergent diagrams in Il37ll starts from the Z2 broken phase and shows that a logarithmic renor- 



maUsation of the bare mass parameter is necessary to cancel divergent contributions, 

Ao = A^ + ^ (In A^ + c) , (57) 

where c may be any constant to fix the renormalisation scale. In analogy to the quenched model this constant is now 
fixed to obtain 

20 



/c 




Figure 16: Extrapolation of the critical couplings of the bosonic 0'* model based on the SLAC derivative to the continuum limit according to the 
functional form (52). 



Figure 17: The divergent Feynman diagrams for the full A^ = 1 Wess-Zumino model in the Z2 symmetric phase. 



with Af^2 defined for the SLAC derivative in Eq. (|54] |. In contrast to the quenched model only a factor '1' in front 
of the divergent part is needed that arises from a partial cancellation of the Feynman diagrams shown in Fig. [17] 
Compared to the Af ~ 2 Wess-Zumino model there is no complete cancellation and a divergence remains. Again the 
given renormalisation procedure amounts to a normal ordering of interaction terms with mass parameter /t in the Z2 
symmetric phase. Eventually a renormalised coupling is defined in the continuum limit similarly as in the bosonic 
case according to Eq. (ISOl l, but now with /tg given through Eq. (l58l l. 

3.2.2. The Pfaffian 

The field content of the M = 1 Wess-Zumino model includes Majorana fermions for which the fermionic path integral 

yields a Pfaffian, 



D^D^Z-c 



-sem-'ip'^Mmi' 



I B4> Pf Af[0]c-^»[ 



(59) 



where M[(j)] is an antisvmmetric matrix! I The Pfaffian is (up to a sign) the square root of the determinant, (Pf Af)^ 
det M, and follows from Grassmann integration, 

1 ^ 



(60) 



CTeS'2 



In practice the Pfaffian is computed as described in ISTllSSll with complexity 0{N^). Similar to the case of Dirac 
fermions in Sec. 12.2. ll this Pfaffian may have a fluctuating sign and Monte-Carlo simulations are carried out with the 
effective action 

5eff = S-B - In |Pf A/| = 5b - i In dct il/ => Z = [ D(j) c^^"'!"^! . (61) 



^'in general, as it is the case e.g. for supersymmetric Yang-Mills theories 186 
antisymmetric but does not need to be anti-Hermitian. 
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I, the matrix 1\I is not necessarily real. However, it is still 



Nevertheless, the sign of the Pfaffian must be taken into account by reweighting of measurements. 
With the chosen representation of the CHfford algebra the fermion matrix is given by 



nSLAC qSLAC _ pi 
^SLAC ^ p/ "^gS 



M = ( .SL^C , w ° ^SLAC ) ■ (62) 



By using general identities for the Pfaffian of an antisymmetric matrix A e ]R2nx2n ^j^^ general matrix B E R2nx2n^ 

Pi{A) = (-)" Pf (A"^), Fi{BAB'^) = dct(B) Pf (A), (63) 

transformation properties under the Z2 symmetry (j> — > —0 of the bosonic potential P^/2 can be derived. For any 
lattice-derivative with antisymmetric matrix representation (e.g. for the SLAC derivative) on a lattice with N points 

di do — P'\ transpose ^ \N r,c f ^^1 ^^0 + P' 






-A (64) 



holds true. Under the Z2 transformation P' — )■ —P' the Pfaffian is invariant on even lattice volumes and hence it de- 
fines a discrete symmetry of the full model. On the other hand, for an odd number of lattice sites the Pfaffian changes 
sign under P' -^ —P'. The seemingly inconsistent behaviour where symmetry properties depend on the number of 
lattice points can be resolved by a closer look at the SLAC derivative. For this derivative the number of lattice points 
is directly related to the boundary conditions of the fields. By imposing the natural condition that the spectrum of the 
lattice derivative operator lies symmetric around the real axis in momentum space an even number of lattice points 
must be used for antiperiodic boundary conditions and an odd number for periodic oneso In consequence the Pfaffian 
changes sign under P' -^ —P' for periodic boundary conditions. Then the contribution of every field configuration 
with positive Pfaffian is cancelled in the path integral by the contribution of the transformed configuration with nega- 
tive Pfaffian and the same bosonic action. This implies a vanishing Witten index, which is nothing but the path integral 
with periodic fermionic boundary conditions. For antiperiodic fermions in the temporal direction the Pfaffian keeps 
its sign under the Z2 symmetry in accordance with the positive definite partition function for the thermal ensemble. 

For Wilson fermions the situation is different. In that case the derivative has still an antisymmetric matrix repre- 
sentation but the contribution of the Wilson term leads to 

di do — P' — §^\ transp. , jv 



-l\ (65) 

10) 

= ^ (-)^Pf 

Changing the sign of P' can in general only preserve the modulus of the Pfaffian if the Wilson parameter r changes its 
sign, too. Since the lattice theory is defined with fixed r the Z2 symmetry will be broken by the Wilson term, similar 
to the M = 2 Wess-Zumino model 16011. For that reason the discretisation based on the SLAC derivative extends the 



-di 
- P' - §A 




ao + p' + §A\ 
d, ) 


do-P'-- 


;a 


do + P' + ^A 
-9i 



works 11371 l80l 18911 by implementing the Z2 symmetry of the continuum model explicitly on the lattice. 



3.2.3. Symmetries, boundary conditions, and ground states 

The static properties under the Z2 symmetry that depend on the boundary conditions for the dynamical fermions can 
be directly related to the ground state structure of the full model. In the case of broken Z2 symmetry one ground 
state is correlated to positive (p while the other one is correlated to negative </>. The analysis of sign Pf M reveals that 
in exactly one of these states the sign flips when one changes the boundary conditions. A change of the fermionic 
boundary conditions amounts to an insertion of (—1)^'' into the path integral which exactly implies the existence of 
one bosonic and one fermionic ground state. 



^"^Nevertheless, the squared SLAC derivative can be defined unambiguously for periodic fields with an even number of lattice points. 
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p(4') 




Figure 18: Probability density of the volume averaged bosonic field for periodic boundary conditions on a 9 X 9 lattice (left panel) and thermal 
boundary conditions on an 8 X 9 lattice (right panel) at coupling / = 100 and A = 0.1. The histograms are computed (with a statistics of 6 ■ 10® 
configurations) separately for fixed signPf M. 

These relations have been checked in the dynamical ensemble with (small scale) lattice simulations with periodic 
and antiperiodic temporal boundary conditions for the fermions at couplings A = 0.1 and / = 100 (see Fig. [TsT l. 
The different boundary conditions at finite temperature imply that the functional forms of both histograms need 
not necessarily coincide. However, all configurations with > keep sign Pf 1\I under a change of fermionic 
boundary conditions, such that the bosonic ground state is located at ^ > 0, while the fermionic has a support < 0. 
These relations have been checked at further couplings / € [10, 100] and sigii0 • sigiiPf A/ > is found on every 
configuration. Similar to the broken supersymmetric quantum mechanics in Sec. |2]it is necessary to use thermal 
boundary conditions to have a positive measurel^j 

This choice of thermal boundary conditions has further implications for the analysis of supersymmetry breaking. 
Apart from the explicit supersymmetry breaking introduced in the lattice theory by the finite lattice spacing and finite 
volume there is a further explicit breaking introduced by the finite temperature. For that reason an analysis of the spoyi- 
taneous supersymmetry breaking in the continuum theory will involve the limit of infinite volume ("thermodynamic 
limit"), vanishing temperature, and vanishing lattice spacing of the lattice theory. 

3.2.4. Z2 breaking 

The numerical survey of the phase diagram starts in analogy to the quenched case with the determination of the phase 
transition associated to the Z2 breaking. The binder cumulant U is thus computed for different volumes at fixed 
lattice spacing with thermal boundary conditions. However, a comparison to results obtained with periodic boundary 
conditions (without reweighting) is enlightening. For A = 0.3 the intersection point of the Binder cumulants is 
independent of the chosen boundary conditions (see Fig.ll9b. This behaviour is explained by the ground state structure 
at infinite lattice volume. In the Z2 broken phase (for small fi^) the system resides in only one ground state with fixed 
signPf M. In that case expectation values are insensitive to the boundary conditions for the fermions and periodic 
boundary may be imposed. For that reason it is quite safe to approach the phase transition from the Z2 broken region 
and extract a critical coupling from a crossing of Binder cumulants at the edge of the Z2 broken phase. But in the 
Z2 symmetric phase reweighted expectation values are undefined for periodic boundary conditions and un-reweighted 
values may only be used with caution. 

The phase transition has also been determined for two further lattice spacings with thermal boundary conditions 
(see Fig. |20] | and the results for every A are in full agreement with a critical coupling of /^ = 21.1(1.1). At this 
point the numerical precision is just not sufficient to resolve any running of the critical coupling with varying lattice 
spacing, mainly because of the accessible lattice sizes. Hence the determined critical couplings are used as continuum 
critical coupling for the Z2 breaking with a broken phase for / > /c and a symmetric phase for f < fc- 



^^The Pfaffian is found to be strictly positive for finite temperature simulations with several couplings on lattices up to 20 X 21. However, for 
even larger lattices the calculation of the Pfaffian is getting extremely time consuming and a representative statistics has not been generated. 
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Figure 19: Binder cumulants for tliermal (left panel) and supersymmetry preserving (right panel) boundary conditions determined at fixed A = 0.3 
and varying fj.^ (with a statistic of lO** configurations). The shaded area denotes the error bounds of the critical p?. 
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Figure 20: Binder cumulants for thermal boundary conditions determined at fixed A : 
indicated by the shaded area. 



: 0.2 (left panel) and A = 0.4 (right panel) with a critical (1^ 



3.2.5. Supersymmetry breaking 

As found in Sec. 13.2.31 the M = \ Wess-Zumino model with the chosen prepotential possesses one bosonic and 
one fermionic ground state which are related through the Z2 symmetry. In the Z2 broken phase one definite ground 
state is selected and this ground state will be the supersymmetric one because its partner state is not present in the 
physical spectrum at infinite volume. On the other hand, in the Z2 symmetric phase there may exist a supersymmetric 
ground state - in contrast to supersymmetric quantum mechanics does the Z2 symmetry not necessarily exclude a 
supersymmetric ground state. Thus it is necessary to study supersymmetry breaking on its own. 

The direct way to study supersymmetry breaking is given by Ward identities that are related to the supersymmetry 
transformation. If there is one broken Ward identity then supersymmetry is broken. The simplest Ward identity which 
is inherently related to the ground state energy is constructed from the transformation of the fermionic field. 



V- 



d X Sifj 



y- 



A'-xP 







with 



<F> 



{V)=G 



{v) = {v)/\/X, 



V 




(66) 



It follows that the dimensionless prepotential V serves as an indicator for supersymmetry breaking. However, super- 
symmetry is explicitly broken by the finite lattice spacing, finite temperature and finite volume. Therefore dynamical 
supersymmetry breaking should be examined in the limit of infinite lattice volume N — > cx), in combination with the 
continuum limit a — > 0. 
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Figure 21: Prepotential ("P) at fixed A = 0.1 and Nf, = 35 for the Z2 broken phase with / = 100 (left panel) and the Z2 symmetric phase with 
/ = 10 (right panel). 
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Figure 22: Extrapolation of (P) to infinite volume at fixed A and couplings / = 100 (left panel) and / = 10 (light panel). Depending on the 
lattice size from 10** up to 5 • 10^ configurations have been used. 

The first limit to be studied is the limit of vanishing temperature at fixed lattice spacing (A — 0.1) and fixed 
spatial volume (iVs = 35). Simulations have been performed at / = 100 (Z2 unbroken) and / = 10 (Z2 broken) and 
the results plotted in Fig.|2T]indicate that for N^ > N^ the finite temperature corrections become negligible in both 
phases. Hence in the following nearly quadratic lattices are chosen with Ni ^ N^ ± lO But already at this point the 
prepotential is one order of magnitude larger in the Z2 symmetric phase than in the Z2 broken phase. 

Now that the finite temperature effects are under control the infinite volume limit can be performed. Here, the 
infinite volume limit is taken before the continuum limit to finally work out the effect of the finite lattice spacing. To 
accomplish an extrapolation the data at fixed finite lattice spacing are fitted with functions of the form 



{V){N,) = A + BN-^ + CN~ 



(67) 



and extrapolated to A^s — > 00. For most of the couplings / lattices with N^ G {25, 27, 31, 35, 43, 63} are used. The 
two examples with A = 0.1 and / € {10, 100} in Fig.l22lillustrate the validity of the chosen extrapolation formula. 

The last limit to be taken is the continuum limit. In the simplest case corrections are of O (a) and a linear extrapola- 
tion to the continuum limit is possible. The extrapolation is done at each coupling / G {10, 12.5, 16, 20, 25, 40, 100} 
and is shown for the extreme cases of the present study, at / — 10 and / — 100, in Fig. |23] where the validity of a lin- 
ear extrapolation is visible. In these cases a complete coincidence between restored Z2 symmetry and spontaneously 
broken supersymmetry is visible. All continuum extrapolated Ward identities in the considered coupling range are 
plotted in Fig.l24land listed in Tab.|3] 



^^The choice of the sign in TVs it 1 depends on the efficiency of the involved Fourier transformations. 
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Figure 23: For the Z2 broken case (/ = 100, left panel) the Ward identity is fulfilled with ("P) = —0.0012(19). For restored Z2 symmetry at 
/ = 10 (right panel) {V) = 0.0444(10) is obtained. 
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Figure 24: Dimensionless prepotential "P over inverse coupling / ^. The region shaded in red indicates the Z2 restored coupling range where 
supersymmetry is broken while the region shaded in gray gives the error bound for the Z2 phase transition. 



The calculation of the dimensionless prepotential that serves as Ward identity by Eq. ( 1661 ) has shown that after tak- 
ing all necessary limiting procedures supersymmetry is broken whenever the Z2 symmetry is restored. Nevertheless, 
one inconsistency shows up for / slightly above the critical coupling for Z2 breaking (approximately / e [/c, 27]). 
In that region the Z2 symmetry is broken while the Ward identity is not fulfilled. This is in contradiction with the 
fact that a broken Z2 symmetry strictly implies a restored supersymmetry after all limits have been taken. Thus an 
analysis of possible systematic errors is in order. 

Firstly the used extrapolation formulae for the specific limits may not be sufficient in every case. Secondly close 
to the critical coupling at fixed A where a second order phase transition (related to an infinite continuum A) with 
diverging correlation length occurs, the considered lattice volumes may still be too small to be in the applicability 
range for an infinite volume extrapolation with Eq. ( l67b . Thirdly the ordering of limits may be of importance. As it 
has been found in the supersymmetric quantum mechanics and the Af = 2 Wess-Zumino model the sign problem that 
arises from the fermionic part of the action becomes worse for larger volumes at fixed lattice spacing and is weakened 
in the continuum limit at fixed physical volume. Although the sign problem is completely absent for small lattices 
signPf M has not been computed for the larger lattices due to the numerical complexity of 0{N^). Therefore sign 
problems cannot be excluded for large lattice volumes and an extrapolation to infinite volume prior to the continuum 
limit may turn out to be insufficient if no reweighting with sign Pf M is done. 
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/ 


{V) 


100 
40 
25 
20 
16 
12.5 
10 


-0.0012(19) 
-0.0039(36) 
0.0032(11) 
0.0093(33) 
0.0150(28) 
0.0277(27) 
0.0444(10) 



Table 3: Dimensionless prepotential V after extrapolation to the limit of vanishing temperature, infinite volume, and vanishing lattice spacing (in 
given order). 







TV, = 49 


N,^63 


N, = 81 


cont. 


LVX = 19.9 


bos. 


0.379(2) 


0.379(4) 


0.383(4) 


0.389(10) 


i\/A = 28.2 


bos. 


0.295(4) 


0.293(5) 


0.293(5) 


0.289(12) 


L\/A = 34.5 


bos. 


0.267(6) 


0.262(8) 


0.261(7) 


0.251(20) 


lVX = 39.8 


bos. 


0.238(10) 


0.228(12) 


0.231(10) 


0.218(30) 


L\/A = 19.9 


ferm. 


0.300(1) 


0.294(2) 


0.291(2) 


0.277(4) 


L^/X = 28.2 


ferm. 


0.276(1) 


0.263(2) 


0.262(2) 


0.237(5) 


lVX = 34.5 


ferm. 


0.270(1) 


0.261(2) 


0.252(2) 


0.225(5) 


Ly/X = 39.8 


ferm. 


0.270(1) 


0.254(3) 


0.245(2) 


0.204(6) 



Table 4: Bosonic (upper rows) and fermionic (lower rows) masses rrtg ;p/ vA for fixed physical volumes LvA and varying lattice spacing together 
with the continuum extrapolation at coupling / = 100. 

Even with taking these possible systematic errors into account a supersymmetry breaking phase transition is 
confirmed where the corresponding critical coupling coincides with that of the Z2 phase transition. Nevertheless, 
the Binder cumulant technique for the Z2 symmetry breaking provides a more reliable way to determine the critical 
coupling because the extrapolation does not directly involve questionable extrapolation formulae. 



3.2.6. Masses 

Ward identities are indicators for the restoration of supersymmetry. Now, after the phase structure of the theory is 
settled, further physical observables are of interest. Amongst them are the particle masses or the energy difference 
between the ground state and the first excited state. Here, one expects a fundamentally different behaviour of the 
masses in the different phases 17711 . In the supersymmetric phase a degeneracy between the (finite) bosonic and 
fermionic mass is expected, similar to the M ~ 2 Wess-Zumino model. For broken supersymmetry a goldstino should 
arise as massless fermionic mode while the physical spacetime volume serves as a regulator for the bosonic mass, 
which itself eventually vanishes in the infinite volume limit. 

The analysis starts with the Z2 broken phase at a fixed coupling of / = 100. In this phase it is necessary to 
project the (finite volume) lattice simulations onto one ground state to mimic the suppression of tunnelling events 
in the infinite volume liniitPJ This is the same technique as used in Sec. |2] for the SQM as well as for the M ~ 2 
Wess-Zumino model in 16011 . and it is necessary to finally extrapolate the obtained masses to the infinite volume 
limit. Since a projection to one ground state is performed it is not required to stick to thermal boundary conditions 
as discussed in Sec. 13.2.41 Thus, in order to remove the supersymmetry breaking introduced by a finite temperature, 
periodic boundary conditions are used also for the fermionic field to obtain correlators and masses in the Z2 broken 
phase. Furthermore only square lattices are investigated to simplify the comparison of different lattice spacings and 
physical volumes. 



Here the configurations are projected without loss of generality to the bosonic ground state. 
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Figure 25: Connected bosonic (left panel) and fermionic (right panel) correlator in the 
LvA = 39.8 for different lattice spacings. 
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Figure 26: Left panel: Masses at finite lattice spacing at coupling / = 100. Right panel: Continuum extrapolated masses together with the error 
bounds of an extrapolation of the fermionic mass to the infinite volume limit (shaded area). 10® configurations have been used for each data point. 



Masses are extracted from the correlators 

CB(t) =iVs"^^ (0(4,^)0(0,,;')) ^"d Cpit) = N-"^ ^ {i'a,{t,x)'^a,{0,x')) 



(68) 



using a cosh fit in a range t E [L/3, 2L/3]. It is obvious that the correlators at fixed physical volume depend on 
the lattice spacing (in Fig.|25]the fermionic correlator shows larger discretisation errors) and extracted masses must 
therefore be extrapolated to the continuum limit. The continuum value is reached via a linear extrapolation that has 
already been used successfully for the continuum extrapolation of results based on the SLAC derivative in unbroken 
supersymmetric quantum mechanics in i42r . The results obtained on lattices with N^ € {49, 63, 81} for the bosonic 
and fermionic masses at four different physical volumes, together with the continuum extrapolation, are given in 
Tab.|4]and are plotted in Fig. |26] (left panel) for finite lattice spacings. Finally an infinite volume extrapolation of the 
continuum results is necessary. Although boson masses approach the fermion masses at larger volumes, as predicted 
by supersymmetry, the statistical accuracy is not sufficient for a reliable extrapolation. Therefore only the fermionic 
masses are extrapolated linearly (see Fig.|26l right panel) to the infinite volume limit, resulting in m.p/y/X = 0.14(1). 
Unfortunately, the bosonic mass can only be assumed to take the same limit. At this point a note concerning the 
finiteness of the bosonic mass is in order. The breaking of the discrete Z2 symmetry does not lead to Goldstone 
bosons and a finite mass is not excluded. On the other hand, in the Z2 symmetric phase the continuous supersymmetry 
will be broken and massless goldstinos should appear in the spectrum of the theory. 

The analysis of masses is therefore continued in the phase with broken supersymmetry and restored Z2 symmetry. 
In this phase it is mandatory to use thermal boundary conditions since bosonic and fermionic ground state participate 
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Figure 27: Fermionic correlators at coupling / = 10 in the thermal ensemble for different physical volumes on a 64 X 63 lattice (left panel) and 
for different lattice sizes on a fixed physical volume Lyk = 44.5 (right panel). For symmetry reasons only the t/L < 0.5 range is shown. 



equally well in the path integral with unsuppressed tunnelling even at infinite volume, 
broken supersymmetric quantum mechanics, show up as massless modes in 



cm = N- 



E 



{i^l,it,x)i^l,{0,x') - '^2,{t,x)ll^2,(0.x')) 



Goldstinos will, similar to 



(69) 



the spinor component combination that yields a cosh form for thermal boundary conditions. At / = 10 the correlator 
has been computed at varying physical volume with fixed lattice size and for fixed physical volume with varying lattice 
spacing (see Fig.|27]i|3 A constant part of the correlator is clearly visible and independent of the physical volume or 
lattice spacing. This is an unambiguous sign of a goldstino. 

To get a complete picture the bosonic masses are calculated as well. These are expected to vanish in the infinite 
volume limit as predicted by renormalisation group methods in 17711 . Again, the connected bosonic correlator at 
/ = 10 is computed and it is found to be composed of a part with nearly vanishingmass, corresponding to the first 
excited state, and a part that arises from higher excited states (see Fig.|28] left panel)|3 The masses of the first excited 
state are now extrapolated (linearly) to infinite volume (see Fig. |28] right panel). Here, only constant lattice sizes 
TV = 64 X 63 are used and the continuum limit is not performed. However, it has been checked for LvA = 44.5 on 
lattice sizes iV = 80 x 81 and A^ = 108 x 105 that the discretisation errors are still below the statistical errors. The 
infinite volume extrapolation is in agreement with m^ ex L^^ with an extrapolated value of me/VA = —0.002(10), 
i.e. the bosonic mass vanishes after the infrared regulator is removed, in agreement with the results obtained by a 
functional renormalisation group approach ll77ll . 



4. Conclusions 

For the case of a supersymmetric quantum mechanics with dynamically broken supersymmetry observables that are 
computed using a lattice regularisation with the SLAC derivative completely coincide with results obtained from the 
diagonalised Hamiltonian. The (bosonic/fermionic) nature of both ground states can be explained with the impact of 
a change in boundary conditions on the fermionic determinant. Correlators computed in the thermal ensemble show 
a constant part for large t, which is a remnant of the degenerate ground states. With a projection to one ground state 
the constant part is still visible in the fermionic correlator, which goes at hand with the massless fermionic excitation 
implied by the degeneracy. In the bosonic correlator the constant part vanishes and the remaining exponential fall-off 
corresponds to the first excited state in the bosonic spectrum. 



^**The non-statistical fluctuations showing up in the con'elator can be traced back to the non-locality of the SLAC derivative. They will decrease 
in the continuum limit, as visible in Fig. l27l (right panel). 

^'The exponential decay of higher excited states is also visible for t\/X < 10 in the fermionic con'elator, cf. Fig. 1271 
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Figure 28: Left panel: Bosonic correlator at couping / = 10 on a 64 X 63 lattice for fixed physical volume LvA = 63 where the contributions of 
first and higher excited states have been separated. Right panel: Bosonic mass of the first excited state on 64 X 63 lattices and the infinite volume 
extrapolation (shaded area) for / = 10. 



On the level of Ward identities it is checked that the ground state is not invariant under the supersymmetry and 
a simple Ward identity is not fulfilled in the limit of vanishing lattice spacing and temperature, as predicted from 
diagonalising the Hamiltonian. In contrast to purely bosonic scalar models, the constraint effective potential does 
not flatten out towards the conventional effective potential, which is explained by the fact that no interpolating states 
between bosonic and fermionic sector are accessibly in the finite temperature path integral. 

We showed that it is possible to analyse supersymmetric quantum mechanics with broken supersymmetry with 
lattice methods based on the SLAC derivative. The physical properties can be determined reliably and accurately. 
Supersymmetric quantum mechanics has been very useful to test sophisticated techniques that find applications in 
higher dimensional models. 

The analysis of the J\f ~ (1,1) Wess-Zumino model in two dimensions aimed at observing and understanding 
dynamical supersymmetry breaking from first principles. A lattice regularisation based on the SLAC derivative is 
used and the choice of this regularisation is justified utilising the quenched model where a complete agreement of the 
obtained critical coupling with the reference value Ii83il is found. 

With this discretisation the Z2 symmetry breaking is analysed and a renormalised continuum coupling is defined. 
For the first time a regulator independent critical coupling is determined from lattice simulations. From the computa- 
tion of a Ward identity a complete coincidence between the restoration of Z2 symmetry and the dynamical breaking 
of supersymmetry is obtained. The computation of masses in the continuum limit for different physical volumes com- 
pletes the analysis and agrees with the picture of a finite and equal bosonic and fermionic mass in the supersymmetric 
phase and the occurrence of a massless goldstino for broken supersymmetry. 

In future works it may be checked on the Ward identities by taking the infinite volume limit after the continuum 
limit has been carried out to suppress possible systematic errors arising from the sign problem. Masses of higher 
excited states could be within reach by using improved correlators. Finallya completely independent calculation with 
different discretisation is desirable. A formulation with Wilson fermions 13711 is a natural choice, in particular because 
efficient lattice methods are already available 19011 . However, in that case one must ensure that a spontaneous breaking 
of the Z2 symmetry is not influenced by the unavoidable explicit Z2 symmetry breaking induced by the Wilson mass 
term. 

In both models the SLAC derivative has proven to be successfully applicable, which can be traced back to the 
absence of gauge fields in these models. Therefore it may be advantageous to consider the SLAC derivative as an 
interesting alternative to Ginsparg- Wilson fermions also in simulations of the four dimensional Af = 1 Wess-Zumino 
model jgifl . 

Supersymmetrically improved lattice actions inevitably include periodic boundary conditions for fermionic fields 
which lead to a severe sign problem in supersymmetric theories with a dynamically broken supersymmetry. There- 
fore the applicability of the improvement programme, where one keeps part of the supersymmetry intact, becomes 
questionable for these theories. It will be necessary to analyse if improved actions exist that give rise to a correct 
continuum limit without fine tuning even for thermal boundary conditions. 
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